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Extensive investigations show that QED 3 exhibits dynamical fermion mass generation at zero 
temperature when the fermion flavor N is sufficiently small. However, it seems difficult to extend 
the theoretical analysis to finite temperature. We study this problem by means of Dyson-Schwinger 
equation approach after considering the effect of finite temperature or disorder-induced fermion 
damping. Under the widely used instantaneous approximation, the dynamical mass displays an 
infrared divergence in both cases. We then adopt a new approximation that includes an energy- 
dependent gauge boson propagator and obtain results for dynamical fermion mass that do not 
contain infrared divergence. The validity of the new approximation is examined by comparing to 
the well-established results obtained at zero temperature. 
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I. INTRODUCTION 

It is well established that (3-|-l)-dimensional quantum 
electrodynamics (QED) can describe the electromagnetic 
interaction between charged elementary particles with 
very high precision. QED dehned on (2-|-l)-dimensional 
space-time, dubbed QED 3 , is safer in the high-energy 
region than its (3-|-l)-dimensional counterpart. In par¬ 
ticular, QED 3 is a superrenormalizable held theory and 
therefore free of ultraviolet divergence. Extensive inves¬ 
tigations have found that QED 3 exhibits a number of 
interesting physical properties such as dynamical chiral 
Symmetry breaking (DCSB) lll- llii. asymptotic freedom 

, and permanent conhnement [III l20l . |2I| . Apparently, 
QED 3 is more like four-dimensional QCD than QED 4 . 
For this reason, QED 3 is often considered as a toy model 
of QCD 4 in the context of particle physics. More in¬ 
terestingly, QED 3 has proven in the past twenty years 
to be an effective low-energy model for several strongly 
correlated condensed-matter systems, including high-Tc 
cuprate su per conductors I22M33I . spin-1/2 Kagome spin 
liquid 134 1351 , graphene [36l438l| , quantum critical sys¬ 
tems [^, and Kane-Mele model with weak-extended 
Hubbard interaction (40l| . 

Although QED 3 can be well controlled in the ultravi¬ 
olet regime, it encounters infrared problems since both 
Dirac fermions and U(l) gauge held are massless. For 
instance, the theory exhibits infrared divergence if the 
free gauge boson propagator is utilized in the perturba¬ 
tive calculations. Appelquist et al. showed that such 
an infrared divergence can be naturally erased by includ¬ 
ing dynamical screening effect of massless fermions into 
the effective gauge boson prmagator 0. Based on this 
scheme, Appelquist et al. Q investigated the Dyson- 
Schwinger equation (DSE) of fermion self-energy func¬ 


tion, and found that the massless fermions acquire a hnite 
dynamical mass, which induces DCSB, when their havor 
is below certain threshold, N < Nc- Most of the existing 
analytical and numerical calculations agree 

that the critical havor is Nc ~ 3.5 at zero temperature. 

In the application of QED 3 to condensed matter sys¬ 
tems, DCSB is usually interpreted as the formation of 
Heisenberg quantum antiferromagnetism [H, Ihol - I^ . It 
was found that DCSB in QED 3 with hnite gauge boson 
mass can still appear if the mass of the gauge boson is 
not large enough . This model can describe the coex¬ 
istence of antiferromagnetism and superconductivity in 
high-temperature cuprate superconductors [4l|. 

The fate of dynamical mass generation at Hnite tem¬ 
perature has also attracted considerable interest. 

Notice that hnite temperature QED 3 is not only inter¬ 
esting from the viewpoint of quantum held theory, but 
of practical value since QED 3 has wide applications in 
condensed matter physics. An important problem is to 
estimate the critical temperature Tc at which dynamical 
mass generation is destroyed by thermal huctuation. 

Unfortunately, the study of dynamical mass generation 
at hnite temperature seems to be more difficult than the 
case of zero temperature. A main obstacle is that the 
Lorentz invariance is explicitly broken at hnite temper¬ 
ature. The integration over 3-momenta k = (fco,fci,fc 2 ) 
that appears in the DSE is replaced by an integration 
over 2 -momenta k = (fci, k 2 ) and a summation over imag¬ 
inary frequency fco = uin, where is the Mastubara 
frequency with n being an integer. As a consequence, 
the DSE becomes much more complicated, so it is usu¬ 
ally necessary to make certain approximations to perform 
algebraic calculations. A number of different approxi¬ 
mations have been proposed to study dynamical mass 
generation in hnite temperature QED 3 IdOl - ls^ . 

The most frequently used one is the so-called instanta¬ 
neous approximation, which assumes that the gauge bo¬ 
son propagator is completely independent of energy. 
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In previous works [42I-I4M Isol - f^ , the transverse compo¬ 
nent of gauge boson propagator, Ay , is usually ignored 
on top of the instantaneous approximation, i.e.. 


A^j.(0,q) Aoo(0,q). 


( 2 ) 


Under these approximations, the DSE for fermion mass 
becomes much simpler and can be easily solved. How¬ 
ever, as showed in Ref. [ 1 ^, the transverse component of 
gauge boson propagator Ay, if incorporated in the DSE, 
induces an infrared divergence at finite temperature. To 
obtain convergent results, Lee proposed to simply neglect 
the contribution from Ay (d^ . 

However, neglecting the transverse component of gauge 
boson propagator is actually problematic. An important 
reason for QED3 to exhibit dynamical fermion mass gen¬ 
eration is that the gauge interaction is long-ranged, which 
is guaranteed by the gauge invariance. If the gauge boson 
acquires a finite mass, the dynamical fermion mass gener¬ 
ation will be strongly suppressed [ 4 l| . At finite tempera¬ 
ture, the Lorentz invariance is explicitly broken, then the 
longitudinal and transverse components of gau ge boson 
propagator behave very differently [ 25 I - IT 7 I . IsTHS^ : the 
former becomes short-ranged due to thermal screening, 
whereas the latter remains long-ranged as a consequence 
of gauge invariance. Therefore, the transverse component 
of gauge boson propagator should play more significant 
role than the longitudinal component at finite temper¬ 
atures. This judgement is supported by the extensive 
recent analysis of nontrivial properties induced by gauge 
interaction [ 25 l - [ 27 L . Overall, the transverse com¬ 

ponent of gauge boson propagator needs to be included in 
an appropriate manner in the study of dynamical fermion 
mass generation at finite temperature, which is the mo¬ 
tivation of our work. 

In this paper, we revisit the issue of dynamical mass 
generation in finite temperature QED3. In order to ob¬ 
tain physically meaningful results, we go beyond the 
widely used instantaneous approximation, and employ a 
new approximation that ignores the energy dependence 
of dynamical fermion mass, i.e., m(poi p) m(0, p), but 
maintains both the longitudinal and transverse compo¬ 
nents of gauge boson propagator. Different form the in¬ 
stantaneous approximation, the new approximation does 
not completely neglect the energy dependence of the 
gauge boson. After performing extensive numerical com¬ 
putations, we find that the dynamical fermion mass is 
free of divergence under the new approximation. We also 
examine the validity of the new approximation by com¬ 
paring to the case of zero temperature, and show that it 
is better than the instantaneous approximation. 

As aforementioned, QED3 can serve as an effective 
low-e nergy f ield theory for several condensed matter sys¬ 
tems 22 h 33 I|. In these systems, there are always certain 


amount of disorders, which couple to massless fermions 
and induce a finite damping rate. Disorders are responsi¬ 
ble to the anomalous behaviors of a plenty of observable 

S antities of interacting systems of Dirac fermions (bOl - 
. Therefore, it is also interesting to study dynamical 


mass generation of QED3 in the presence of disorders. 
To describe the influence of disorders, a finite fermion 
dam ping rate T is introduced to the fermion propagator 
Ref. . As will be shown below, T plays analogous role 
as temperature T, so its effect on dynamical mass gener¬ 
ation can be studied in a similar manner to the case of 
finite temperature. 

The rest of paper is organized as follows. In Sec. [Hi 
we give the Lagrangian and the relevant propagators. In 
sec.nni we show that the dynamical mass is divergent 
under the instantaneous approximation if temperature or 
damping rate is finite. In Sec. El we solve the DSE for 
dynamical mass by invoking a new approximation, and 
find that the results are free of infrared divergence. The 
nature of the infrared divergence is also discussed. In 
Sec.El we examine the validity of the new approximation 
by comparing to the case of zero temperature. In Sec. lVIl 
we summarize the main results. Detailed calculations of 
polarization functions are given in Appendix 


II. MODEL AND FEYNMAN RULES 


The Lagrangian density for QED3 with N flavors of 
massless Dirac fermions is given by 


N 


1 


C = + e4) - -F^^. 


( 3 ) 


The electromagnetic tensor is related to vector potential 
as — d^A^^. The fermion is described by 

a four-component spinor field tjjj whose conjugate spinor 
field is t/j = The 4 x 4 gamma matrixes are de¬ 

fined as (70,71,72) = (*0-3,icTi, 10-2) ® 0-3, which satisfy 
the standard Clifford algebra {7^,7;^} = 251^1, with the 
metric being = diag(— 1 , — 1 , — 1 ). In ( 2 -l-D dimen¬ 
sions, there are two chiral matrices 73 and 75 (2.f6^. 


73 = * 


0 I 
-I I 


75 = ^ 


0 I 
I 0 


( 4 ) 


which anticommute with 70, 71 and 72. The Lagrangian 
shown in Eq. @ respects a continuous U (2 A) chiral sym¬ 
metry Ip —^ where 9 is an arbitrary constant. 

Once a hnite fermion mass is dynamically generated, the 
global U( 2 N) chiral symmetry is spontaneously broken 
down to its subgroup U{N) x U{N). In this paper, we 
consider a general large N and perform perturbative ex¬ 
pansion in powers of 1 /A. For convenience, we work in 
units with h = kg = 1 and set the velocity zip’ = 1. 

In the Euclidian space, the free propagator of massless 
fermions is 


Go{k) = 


1 




( 5 ) 


at zero temperature. The free propagator of gauge boson 
can be written as 




( 6 ) 
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in the Landau gauge. After including the dynamical 
screening effect due to fermions, the effective gauge boson 
propagator takes the form 

^ + n(g) 

To the lowest order of 1/A^-expansion, the polarization 
function can be obtained from the polarization tensor 
through the relationship 

/ d^k 

^Tr[Go(A:) 7 ^Go(fc + g) 7 .] 

= n(g) , ( 8 ) 

where n(g) = aq with a = Q. 

Pisarski [I| first carefully studied DCSB in QED 3 by 
means of non-perturbative DSEapproach, and showed 
that DCSB takes place for any finite flavor N. How¬ 
ever, subsequent analysis of Appelquist et. al. Q found 
that DCSB can occur only when the fermion flavor N 
is smaller than a critical value Nc, which is Nc = 32/7r^ 
to the lowest order of 1/A^-expansion. Nash then ex¬ 
amined the effect of next-to-leading order correction and 
obtained a critical flavor Nc = (4/3)32/7r^. Pennington 
et. al. [3, @ included the wave renormalization function 
and claimed that DCSB can be realized for any flavor 
iV, although the corresponding dynamical fermion mass 
decreases exponentially with increasing N. Nevertheless, 
their analysis ignored the influence of wave renormaliza¬ 
tion function on the polarization. Later, Maris Q studied 
DCSB by solving a set of self-consistent DSEs for fermion 
propagator and polarization, which contain vertex cor¬ 
rections and therefore satisfy the Ward-Takahashi (WT) 
identities. The calculations of Ref. @ arrived at a finite 
critical flavor Nc ~ 3.3, which is close to that of Ref. Q. 
The key difference between the treatments of Penning¬ 
ton et al. 0 and Maris Q is that the latter included 
the interaction correction to the polarization whereas the 
former did not. It turns out that an appropriate approxi¬ 
mation plays a crucial role in the DSE analysis of DCSB. 
More refined calculations of Fisher et al. revealed a finite 
critical flavor 4jM . The gauge invariance of Nc is 
also discussed [J, [ill [l3 ■ Apart from the DSE approach, 
this problem can be studied by renormalization group 
method, which found that 3 < Nc < i [I^. The criti¬ 
cal flavor Nc obtained in lattice Monto Carlo simulations 
[ 15 I [ 1 ^ is much smaller than that obtained by means of 
DSEs. However, Gusynin et al. argued that the differ¬ 
ence is attributed to the finite volume effect introduced 
in lattice simulations 0 . It is fairy to say that, although 
there is still some debate [ill US, [l^ > most studies have 
obtained a finite critical flavor for DCSB in QED 3 , which 
is roughly Nc ~ 3.5 at T = 0. 

As we go to finite temperature, the dynamical fermion 
mass is expected to be rapidly suppressed by thermal 
fluctuation. The temperature scale Tc at which the dy¬ 
namical mass vanishes defines the critical temperature. 


At T 7 ^ 0, we write the fermion propagator in the stan¬ 
dard Matsubara formalism as 

G(fco,k) (fcg + rsgn(/ep))-),o + ^ . k -f Too’ 

where fco = (2n-|- l) 7 rr with n being an integer. Here, we 
introduce a constant P to represent the fermion damp¬ 
ing rate generated by disorder scattering. This quantity 
measures the strength of the fermion damping effect. For 
more explanation of the origin and the physical effect of 
the constant P, please see Ref. [i^. To the leading order 
of 1 /A^-expansion, the polarization tensor is given by 


]Ve 2 f 

n;,,.(go,q,'7’,mo,r) = ^ / -^;^Tr[G(fco,k) 

X lp.G{kQ -f go, k -t- q)7^], (10) 

where go = ‘i.n'iTT with n' being an integer and (3 = ^. 
The effective propagator of gauge boson now becomes 


A^i.(go,q) = 


A 


fiiy 


9o + q^ + nA(go,q) 

_i__ 

9o+ q^ + nB(go,q)’ 


where 




fc <7^ f. qoq,^\ 




til I Ulj 


_ qiqj \ 




9 I 

Ap,i, and B^i, are orthogonal and satisfy 


Afj,i/ Bp^i/ — 5 fjiv 


qti.qu 


The functions H^i and Hs are defined as 


IIa — ^Hoo, 

q 2 


Ab = H,, - 4noo. 

q 2 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 


The calculational details of the relevant polarization 
functions are shown in Appendix 

In the following, we will consider two approximations, 
namely the popular instantaneous approximation and a 
new approximation to be explained below. We neglect 
the energy-dependence of the polarization functions and 
also the feedback of dynamical mass to the polarization 
functions. Firstly, in the limit that go = 0, P = 0, and 
mo = 0 , the polarization functions are 


nA(q,r) 

nB(q,r) 


16aT 


dx In 


2 cosh 


y/x{l - x)q^ 
2T 


^ Jo 


X tanh 


yjx(l - a:)q^ 
2T 


(16) 


(17) 
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We have used a = —. Secondly, in the limit that go = 
0, T = 0 and toq = 0, the polarization functions are 


n^(q,r) = ^ rin(^)+r 


1 + 


v/q 2 + 4r2 

2|q| 


X in 


y^q 2 + 4r2 _ |q| \ 

\ A/q2 + 4r2 + |q|/ 
+ / dxy/ x{l — 

Jo 

yjx{l - a;)q2 


X arctan 


r 


(18) 


nB(q,r) = ^ 

■TT-^ 


/ c?a;\/a;(l — a;)q 2 

Jo 


X arctan 


\/a;(l - a;)q 2 


(19) 


From these expressions, we can see that the fermion 
damping rate F plays a very similar role to temperature 
T, which allows us to study the effects of temperature 
and fermion damping using the same scheme. To simplify 
later calculations, we can further approximate Eqs. dm) 
and (fT71) by [ij 


nA(q,r) Ri a 


|q| + ciT exp - 


ciT 


nB(q, T) R! Q;|q|tanh 


C 2 |q| 

T 


( 20 ) 

( 21 ) 


where ci = 161n2/7r and C 2 = 2 / 37 r. Analogously, 


Eqs. (HI and (HU) can be approximated by 


, . 16aF / A 

n^(q,F) R p 


HdS! arctan''’™ 


nB(q,r) arctan''’™ 


( 22 ) 

(23) 


where C 3 = 4 / 37 r. It can be checked numerically that 
expressions Eqs. (ISol-dSSI) are very good approxima¬ 
tions for both the high- and low-momentum behaviors 
of Eqs. (fTO - dTOl) . In the following sections, we will use 
Eqs. (I2n)l - (j23p to analyze and numerically solve the DSEs 
for dynamical fermion mass at finite temperature and/or 
finite fermion damping rate. 


III. DSE UNDER INSTANTANEOUS 
APPROXIMATION 


In this section, we present the DSE for dynamical mass 
under the widely used instantaneous approximation. We 
will show that the solutions of DSE are divergent in the 
infrared region whenever T ^ 0 and therefore ill-defined. 
Such an infrared divergence also exists if T = 0 but F 7 ^ 0. 

The free fermion propagator is 


Go(fco,k) = 


(fco -hrsgn(fco))7o + 7- k' 


(24) 


Due to gauge interaction, the fermion may become mas¬ 
sive and the propagator is renormalized to 


G(fco,k) = 


(fco + Tsgn {ko))^o -I- 7 • k -I- m{ko, k, T, F)' 


(25) 


Here, to the lowest order of 1/A'^-expansion, we neglect the wave renormalization function. Now, Go(fco)k) is related 
to G(fco; k) through the following DSE, 


where 


G \po,p) = Gq^{po,p) +m{po,p), 
d^k 


= ( 2 ,)C 

kpj 


7^GF(fco,k)7i2A ( 90 , q) 


(26) 


(27) 


with qo = Po — ko and q = p — k. Substituting Eqs. (IMl) and (1^ into Eq. (I?7l) . and then giving an explicit T- and 
F-dependence to the fermion mass, we obtain 


d^k 


m{ko,k,T,T) 


iiPo,PXX) (3^7 (27r)2 (fco-f rsgn(fco))^ + k2 H-m2(fco,k,T,r) 


1 


1 


9o+q^+n^(90,q,T,F) gg-Hq2-HnB(go,q,T,F)_ ' 


X 


(28) 
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In previous studies the key assumption beneath the instantaneous approximation is to completely 

ignore the energy-dependence of gauge boson propagator. Apart from this approximation, the transverse component 
of gauge boson propagator is widely neglected, i.e., 

q) ^oo(0, q)- (29) 

Now we follow the instantaneous approximation and drop the energy dependence of g au g e boson propagator. Never¬ 
theless, as illustrated in recent works on the nontrivial properties of QED 3 (25l - [^ ISTHS^ . the transverse component 
of gauge interaction plays more important role than the longitudinal component at finite temperature, since the latter 
becomes short-ranged after acquiring an effective thermal mass proportional to T. If T = 0 and F 7 ^ 0, the longitudi¬ 
nal part of gauge interaction also becomes short-ranged due to static screening caused by disorder scattering. In any 
case, the gauge invariance ensures that the transverse component of gauge interaction is strictly long-ranged. It is 
therefore not appropriate to neglect the transverse component. For completeness, here we keep them both and write 
the gauge boson propagator as 


^t^u{0,q) = 


Sf ti (Sr 


QiQj 


2 I 


q 2 -hn^( 0 ,q) q 2 -H 115 ( 0 , q) 


The corresponding DSE becomes 

o 2 


m 


ko 


f (Pk TO(k, T, F) 

1 1 

/ (27r)2 [ko + Fsgn(fco))^ -b k2 -b m2(k, T, F) 

[q2+nn(0,q,T,F) q2 + 05 ( 0 , q, T, F) J 


(30) 


.(31) 


Once the instantaneous approximation is adopted, the dynamical mass m completely loses its dependence on energy. 
Therefore, the frequency summation in the DSE can be written as 




1 


|^(fco+rsgn(fco))^-f k 2 -f m 2 (k,T,r) 


27r 


2 00 


E 


n— — oo 


1 


(n-I- 5 -I-Xsgn (n-I-^)) +Y'^ 


(32) 


where X = ^F, and Y = + m^(k, T, F). With the help of an identity 


SiX,Y) = J2 


n—0 


(n + A )2 -p y 2 2Yi 


= [V'(^ + - m - ^Y)] 


where 'i/' {z) is Digamma function, we get 


= 


■E 


1 


[(n+i+A)%F 2 

Carrying out frequency summation leads to 

So f d^k TO(k, T, F) 


27r2y 


Im 


^(-+X + tY 


m(p,r,F) = 


Im 


(2^)2 ^k2+m2(k,T,F) ^ 

1 1 




1 


2 2ttT 


yk2-f m2(k,r,F) 
2ttT 


-b 


Lq2 + nn(0,q,r,F) q2 + n5(0,q,T,F)J ■ 

This equation can be further simplified. First, we assume F = 0 and T ^ 0, and then rewrite the DSE as 

m(k, T) 


4a /• (fk 
m{p,T) = — 


N J (27r)2 y^k2 + m2(k, T) 


tanh 


y/k^+m'^{k,T) 


2T 


1 


-b 


Lq2+n^(q,T) q2 + nB(q,r)J’ 


where Im 


^ (2 2 ^) ~ f is used in the derivation. Secondly, at T = 0 and F 7 b 0, we have 


(33) 


(34) 


(35) 


(36) 


m(p,F) 


8 a 

"/V 


(Pk TO(k, F) 1 / -b nPCk, F) 

7 —TTT — , — arctan - 

(27r)2 y^k2 -bm2(k,F) tt F 


1 

q2 -bnn(q,F) 


1 

q2 +n5(q,F) ’ 


(37) 
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FIG. 1; Relationship between m{fi)/a and fi/a for different temperatures at F = 0 with = 2, 3 in (a) and (b). 




FIG. 2: Relationship between m{fi)/a and fi/a for different F at T = 0 with = 2, 3 in (a) and (b) respectively. 


where lim 7 ’_>.o Im 


'^{h + IFt ) 


arctan The similarity between the above two equations can be clearly seen. 

At T = 0, the constant fermion damping rate F can be considered as certain effective temperature. This is an 
important reason for us to discuss these two cases simultaneously in this paper. 

To solve the DSE, it is convenient to introduce an infrared cutoff /r. This amounts to assume that the system 
has a finite volume [l^. [^. For finite the dynamical fermion mass obtained from DSE is always free of infrared 
divergence. Our aim is to examine whether the dynamical mass is still well defined as /i —>■ 0. We employ an infrared 
cutoff Hi for momentum |k| and replace |q| with |q| + H 2 - Usually, fii and H 2 should satisfy Hi ~ M 2 - For simplicity, 
we assume Mi = M 2 = M- Now the DSE at F = 0 and T ^ 0 becomes 


i(p,T) = 


Ntt^ 


n2-K 

d|k||k| / dip 
Jo 


i(k,T) 


0 ^k^+m^k,T) 


: tanh 


v/k2+m2(k,T) 

2T 


+ 


(|q| + m)^ + n^dql + M; T) (|q| + /i)^ + Fl^dql + Ml T) 


(38) 


where ip is the angle between k and p. As usual, we choose A = The DSE at T = 0 and F ^ 0 can be transformed 
analogously. 

For different finite temperatures, the dynamical mass at the lowest momentum h with N = 2 and N = S are shown 
in Figs. [T] (a) and (b) respectively. As /r —>■ 0, m{H) does not converge to a finite value, but is divergent. As T further 
decreases, m{H) grows more slowly, but still divergent at finite T. For different values of F, the relations between 
m{H) and h with N = 2 and N = 3 are displayed in Fig.j^Ka) and (b) respectively. Numerical calculations show that 
m{H) is ill defined as ^ > whenever F 7 ^ 0. Apparently, the instantaneous approximation leads to divergent results. 
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IV. DYNAMICAL MASS GENERATION 
UNDER A NEW APPROXIMATION 


We have shown in the last section that the widely used 
instantaneous approximation leads to infrared divergence 
in the dynamical fermion mass obtained at finite temper¬ 
ature or in the presence of a finite fermion damping rate, 
when both the longitudinal and transverse components 
of gauge boson propagator are adopted. It is therefore 
necessary to go beyond this approximation and seek a 
new approximation that could yield convergent results. 
We find it very helpful to employ a treatment utilized by 
Ref. in the study of excitonic insulating transition in 
graphene The key assumption of this treatment 

is to neglect the energy dependence of dynamical fermion 
mass, which was shown [67l| to generate physically reli¬ 
able dynamical mass for Dirac fermions. 

Let us start from the DSE given by (28). Now we 
assume the dynamical fermion mass does not depend on 
energy, namely 


m(po, p,T, r) -)> m(p, T, r). (39) 


After making this approximation, the DSE (1^51) can be 
written in the following form 


w(p,T,r) = 


8 a 


E 

ko • 




{ko + Esgn (fco))^ -I- -f m2(k, T, L) 
1 


Lfc2 + q2 + nA(fco,q,r,r) 


1 


(40) 


^0 "k q^ + nB(fcoi Q) 2^) r)_ 

Notice the effective gauge boson propagator still retains 


an explicit energy dependence, which is thus different 
from the instantaneous approximation. However, the 
summation over ko cannot be performed exactly due to 
the complicated ko dependence of polarization functions. 
To further simplify the DSE, we drop the energy depen¬ 
dence of polarization functions, i.e., 

nA(A:o,q,T,r)^n^(0,q,T,r), (41) 

nB(fco,q,r,r)^nB(0,q,T,r). (42) 


The DSEs then can be formally written as 


m(p,r) 


8 a f m(k, T) 

(27r)2fc2-t-k2 4-m2(k,T) 

A'O 


1 

k^ + (f +UAiq,T) 


1 

+ fc2 + q2+nB(q,T) 


(43) 


in the case oi T ^ 0 and T = 0, and to 


8a I u,n,o 

*(p,r) = ^ / — 


dko 

27r 


d^k 

m(k, I) 


(ko + Tsgn (fco))2 -k k2 -f m2(k, T) 

1 


ko + + n^(q, r) 

1 

A:g + q2-f nB(q,r) 


(44) 


in the case of T 0 and T = 0. The infrared behaviors 
of these two equations will be further analyzed later. It 
is now straightforward to sum over ko , which leads us to 
the following expressions: 


m(p,T) 


m(p,r) 


4a 
'N , 

X 


d^k 


m(k, T) 


-n^(|q|,T)^k 2 -m 2 (k,T) 
' ^k^ + m^{k,T)\ 




^(k,r) 


: tanh 


2T 


v/q 2 +n^(q,T) 


: tanh 


Vq 2 +nA(q,T) 

2T 


+ [nA(q,T)^nB(q,T)]}, 

8a f d^k f m(k, r) 

^ J (27r)2 j (r2 + k2 + m2(k, T) - q2 - n^(q, T))" + 4r2 (q2 + n^(q, T)) 


^ >-k^-m^(k,r)+q^+n^(q,r) fn . f T 

v/k2+m2(k,r) 1^2 I^v'k2 + m2(k,r) 

, 7r r^ + k^ + m"(k,r)-q^-n^(q,r) /^ r^ + k^ + m"(k,r) V 

^2 ^q2 + n^(q,r) ""I q2 + nA(q,r) )_ 

+ [nA(q,r)->nB(q,r)]}. 


(45) 


(46) 












































FIG. 3: Relation between m{^)/a and /i/a for different T with = 2, 3 in (a) and (b). Results are obtained under the new 
approximation. 



FIG. 4: Relation between m(/i)/a and /i/a for different F with = 2,3 in (a) and (b). Results are obtained under the new 
approximation. 


Following the approach used in Sec. IIIII we introduce 
an infrared cutoff /i to Eqs. (l46l) and (l46l) . and then solve 
these equations numerically. In the case that T 0 and 
r = 0, the dependence of dynamical mass m(/i) on /r at 
finite T with N = 2 and N = Sis shown in Figs. [3] (a) and 
(b) respectively. As the infrared cutoff vanishes, /r —?► 0, 
m{^) is saturated to certain finite values, which means 
the infrared divergence encountered under the instanta¬ 
neous approximation does not exist under the new ap¬ 
proximation. Analogous calculations can be carried out 
in the case that T = 0 and F ^ 0. We display the depen¬ 
dence of m{fj,) on /j, with N = 2 and A^ = 3 in Figs. HJa) 
and (b) respectively. It is easy to see that m(/r) also 
approaches finite values as /r —^ 0. Apparently, the new 
approximation adopted in this section leads to convergent 
results for the dynamical fermion mass, and is therefore 
more reliable than the instantaneous approximation. We 


fix the infrared cutoff at /i = 10“^^ and calculate the dy¬ 
namical gap m(|p|) under the new approximation. The 
calculated m(|p|) under the new approximation at differ¬ 
ent values of T or different values of F is shown in Figs. [5] 
(a) and (b) respectively. 

Let us now analyze the origin of the infrared diver¬ 
gence. At T ^ 0 and F = 0, the DSE has the form 


'm{Po,P,T) 


^ qQ—2n7zT '■ 


(Pq m{ko,'k,T) 

(27r)2 k^+lP + rrP{k[), k, T) 


q^ + q^ + nAiqo,q,T) 


q^ + q^ + UB{qo,q,T) 


(47) 


Since m is finite after dynamical mass generation and 
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11 ^ 1 ( 0 , 0,T) oc T in the low energy limit, the terms ap¬ 
pearing in the first and second lines of the integral ker¬ 
nel are safe in the infrared region. However, the term 
appearing in the third line contains a potential infrared 
divergence. To make this transparent, we divide the sum¬ 
mation over ko as follows. 


h 


E 

qo=2nT 


(f(\ 1 

(27r)2 q2 + Ub (go, q, T) 


E 

< 70 —2nT(n/0) 


(P(\ 1 

(27r)2 g2 q2 + Ub (go, q, T) 


(Pq 1 

(2^)2q2+nB(q,T)- 


(48) 


Notice that go = 2nT{n ^ 0) is always finite at finite 
T whenever n ^ 0 , so the first term does not yield any 
infrared divergence. On the contrary, the second term 
is dangerous because nB(q,T) = with m oc ^ for 
small momenta. We now simply focus on the potential 
divergent term, and find that 

f 1 ^ f 

J (27r)2 q2 + nB(q,T) J (27r)2 cq^ 

cJ^ hi c vm/ 

where c = 1 -I- oi. It is clear that this term is divergent 
as the infrared cutoff p, —?► 0. From the above analysis, 
we see that the infrared divergence of DSE comes from 
the zero-energy transfer processes mediated by the sin¬ 
gular transverse component of gauge boson propagator. 
Lee [i^ noticed the existence of infrared divergence, and 
then simply neglected the transverse component of gauge 
boson propagator. The same strategy is widely utilized 
in other works Isol - lssjl . 

More recently, Lo and Swanson @ also stressed the 
existence of infrared divergence and proposed to remove 
this divergence by choosing an appropriate temperature 
dependent gauge parameter. Their approach is basically 
equivalent to considering the following DSE, 


Under the new approximation, the DSE of dynamical 
mass is simply Eq. (1431) . The denominator of the kernel 
of this equation contains a factor of A:g = (( 2 n -|- l) 7 rr)^, 
whose minimum is We may consider this term 

as an effective thermal mass of gauge boson, i.e., rua oc 
, which serves as an infrared regulator and elimi¬ 
nates the potential infrared divergence of the dynamical 
fermion mass. This thermal mass exists only at finite T, 
and vanishes naturally as T —> 0. To examine to what ex¬ 
tent the new approximation is valid, we will show in the 
next section that Eq. (|43)) leads to results qualitatively 
consistent with those obtained at T = 0 Q. 

We then turn to the case of finite fermion damping 
rate. At T = 0, the DSE is given by Eq. (l44l) . The poten¬ 
tial infrared divergence can only come from the second 
term in the bracket. The most singular part is repre¬ 
sented by 


h 


f ^ f cl^q _1_ 

J 2tt J (27r)2 g2 _^q2+nB(go,q,r)' 


In order to analyze potential infrared divergence, we re¬ 
place nB(go,q, r) with nB(q, T), then 


h 


f ^ [ ci^q _ I _ 

J 2Tr J (27r)2 _,_q2 4.n5(q,r)' 


For small momenta, nB(q, F) behaves as nB(q, F) = 
02 q^ with 02 oc p. It is clear that the potential infrared 
divergence can be represented by 


h 


f ^ 

J 2tt 

I 

47r-\/U 


cPq 


(27r)2 g 2 c'q2 


d|q| = 




1 

27? 

(A-/r), 


(Pq 1 

(2?7M 

(51) 


where c' = 1 -I- 02 . As /r —>• 0, /2 is definitely not diver¬ 
gent. This property should be fulfilled no matter what 
approximation is used to calculate the DSE. 

Under the instantaneous approximation, the DSE is 
represented by Eq. dSZl), which contains the following sin¬ 
gular contribution, 


to(po,P,7’) 


k (Pq m(A:o,k,T) 

13 J (27r)2 fcg + k2-I-m2(fco,k, T) 


h 


[ 1 

J (2^)2q2+nB(q,r) 



X 


E 

qQ—2n'KT 


I 

go + q^ + nA(go,q,T) 


+ E 

qo—2n7rT{n^O) 


1 

g2-tq2-tnB(go,q,T) 


This function is divergent as ^ 0. Apparently, the 

instantaneous approximation brings an artificial infrared 
divergence that should not exist. 

Under the new approximation proposed by us, the DSE 
is given by Eg. (Hi)) . The only possible singular part can 
be simply written as 


which ignores the zero frequency (n = 0 ) contribution 
of the transverse component of gauge boson propagator. 
Careful numerical computation of this equation is inter¬ 
esting, but challenging since it is hard to sum over n and 
at the same time integrate over q with high precision. 
This is subjected to future investigation. 


r ^ ^ [ Eq _ 1 _ Kai 

J 2tt J (27r)2/c2 q2-H nB(q,r) 

^ = . ^y 7 (A-M), (54) 

47 rV c' 

which is not divergent and well consistent with Eq. (51). 
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FIG. 5: Momentum dependence of m(|p|) under the new approximation for various values of (a) T and (b) F with p = 10 




FIG. 6: Dependence of on p/a at T 


0 (a) neglecting and (b) including the feedback of Jn(p) to the polarization. 


V. VALIDITY OF NEW APPROXIMATION 


In Sec.IV, we have adopted a new approximation to 
study the DSE for dynamical fermion mass in the pres¬ 
ence of finite temperature or finite fermion damping. 
Compared to the popular instantaneous approximation, 
the main advantage of the new approximation is that, it 
retains both the longitudinal and transverse components 
of gauge boson propagator, and at the same time leads to 
physically meaningful, convergent results. To further see 
this point, it is now interesting to make a more straight¬ 
forward comparison between these two approximations. 

In this section, we consider zero temperature QED 3 in 
the clean limit, namely T = E = 0. This model has al¬ 
ready been extensively investigated, and the dynamical 
mass obtained from DSE is free of infrared divergence, 
which then can be considered as a reference to examine 
the reliability of the results obtained at finite tempera¬ 
ture. For an approximation to be reliable, it should work 
well at both zero and finite temperatures. As T —>■ 0, 
m{T) oughts to approach a well-defined quantity rn{0), 
which should be free of infrared divergence and as close 


in quantity as possible to that obtained directly at T = 0 . 
We now examine whether the new approximation is valid 
according to this criterion. 

If one neglects the feedback of fermion mass to the 
polarization functions, the DSE to the lowest order of 
1/V-expansion is known to be Q 


, . 4a f (Pk m{k) 1 

= :V y k^ + m?{k) g2-bn(g)’ 

with n(( 7 ) = aq. This equation was first solved in Ref.Q, 
and the solution is very well known. Under the instanta¬ 
neous approximation, the DSE is simplified to 


to(p) 


8a 

'n 


m(k) 1 

(27r)2 y/k2 -b m(k) + n(q) ’ 


with n(|q|) = a|q|. Under the new approximation, the 
corresponding DSE is 


m(p) 


4a f cP'k m(k) 1 

^ y (27r)2 ^k2 + m2(k) 

1 

\/k2 -1- TO2(k) 4- yjcp + n(q) ’ 


(57) 
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where n(q) also equals to Q;|q|. The dynamical fermion 
mass as a function ^ obtained in three cases are 

shown in Fig. [6] (a), represented by lines with different 
colors. We notice that m{^) is saturated to finite val¬ 
ues in all these three cases when /r —0. However, the 
fermion mass obtained in the new approximation is closer 
to that obtained directly from Eq. (1551) than the instan¬ 
taneous approximation. It seems that both of these two 
approximations lead to convergent results for dynamical 
mass. Nevertheless, we still need to examine whether 
these results are robust against higher order corrections. 

We then include the feedback of dynamical fermion 
mass to the polarization. The polarization appearing in 
Eq. (|55ll becomes 


n(g,mo) = 




’mg - Ami 

2 g 2 


4g3 


X arcsin 




q2 + 4 too 


(58) 


and the polarization appearing in Eqs. (I56p and (1571) can 
be obtained by replacing q of n(g, mg) with |q|. In addi¬ 
tion, Too is substituted by m{qi). The dependence of dy¬ 
namical mass m{qi) on /i is depicted in Fig.| 6 Kb). We find 
that m{fi) diverges as /i —> 0 under the instantaneous ap¬ 
proximation. Nevertheless, m{fi) obtained in the other 
two cases does not exhibit infrared divergence and re¬ 
mains finite as ^ 0. These results further demonstrate 

that the instantaneous approximation yields unphysical 
results, and that the new approximation is more reliable. 


VI. SUMMARY AND DISCUSSION 

In this paper, we have studied dynamical fermion mass 
generation in QED 3 after including the effects of finite 
temperature or finite fermion damping rate. Many previ¬ 
ous DSE calculations of dynamical fermion mass adopted 
an instantaneous approximation, which is often accom¬ 
panied by simply ignoring the transverse component of 
gauge interaction |42h451 [SOj-lSSj . As already explained in 
the context, at finite temperature or at a finite fermion 
damping rate, the longitudinal component of gauge in¬ 
teraction becomes short-ranged due to static screening, 
whereas the transverse component of gauge interaction 
remains long-ranged as required by the local gauge invari¬ 
ance. It is therefore not appropriate to ignore the more 
important contribution of gauge interaction. However, 


we have showed that, if one adopts the instantaneous ap¬ 
proximation and meanwhile includes the complete gauge 
boson propagator, the dynamical fermion mass exhibits 
infrared divergence. 

We have revisited this problem and employed a new ap¬ 
proximation to calculate the DSE for dynamical fermion 
mass. Under the new approximation, both the longitudi¬ 
nal and transverse components of gauge interaction are 
incorporated, and the results obtained from in the DSE 
are free of infrared divergence. To further examine the 
validity of the new approximation, we have also make a 
comparison to the results obtained directly at zero tem¬ 
perature. In summary, our calculations have showed that 
the new approximation leads to more reliable results for 
the dynamical fermion mass than the widely used instan¬ 
taneous approximation. 

The existence of infrared divergence is not special to 
the issue of dynamical mass generation in finite temper¬ 
ature QED 3 . Indeed, similar divergence appears in a 
number of interacting gauge field theories. For instance, 
Lee calculated the fermion damping rate caused by gauge 
interaction within an effective non-relativistic U(l) gauge 
field theory, and found non-Fermi liquid behavior at zero 
temperature Nevertheless, the fermion damping 

rate diverges at finite temperature [t^. Recently, anal¬ 
ogous divergence is also found in QED 3 defined at finite 
temperature and finite chemical potential [s^. We hope 
the approach proposed and used in this paper could pro¬ 
vide useful insight into this problem. 

Confinement is an important feature of QED 3 . It is 
known from previous studies [n|,[23,[2H that whether this 
model is confining depends crucially on the behavior of 
the polarization function in the low energy regime, which 
is in turn determined by the dynamical fermion mass. It 
would be interesting to apply the new approximation pro¬ 
posed here to carefully calculate the polarization function 
at finite T by including the impact of dynamical fermion 
mass, and then to evaluate the critical temperature for 
the confinement-deconfinement transition, following the 
schemes presented in Refs. (Til. Sim. It would also be 
interesting to examine whether confinement and dynami¬ 
cal fermion mass generation take place simultaneously by 
analyzing the behavior of wave function renormalization 
0 . To address these issues, one needs to incorporate the 
wave function renormalization and the vertex functions 
in the DSEs, which are subjected to future research. 

The authors acknowledge the financial support by the 
National Natural Science Foundation of China under 
grants 11274286, 11174290, and U1232142. 


Appendix A: Calculations of Polarization functions 


In the Appendix, we present the detailed calculations for the polarization functions in the presence of finite temper¬ 
ature T and finite fermion damping rate F. The calculations are performed within the standard Matsubara formalism 
for finite temperature quantum field theory. 
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1. Expression for general ko and general F 


Starting from the effective fermion propagator given by Eq. we write the polarization functions IIoo and 11 ^^ in 
the following forms 


noo(go,q,7’,?7ro,r) = 


Ne 


2 +o“ 

E 

^ n= —oo ‘ 

ni,(qo,q,7’,mo,r) = ^ 

^ n— — oo ‘ 


w? 


Tr [G(fco, k)7oG(fco + 90 , k + q) 7 o] 
Tr [G(fco, k)7iG(fco + qo,k + q) 7 i] , 


(Al) 

(A2) 


where /cq = (2n + l) 7 r //3 and Qq = 27rnY/3 with n and n' being integers, to the leading order of 1/iV-expansion. 
Substituting Eq. ([5]) into Eqs. (ED and (IA2I) . and then using the Feynman parametrization formula 


= dx 


AB Jo [a;A + (1 - x)B]^ ’ 


(A3) 


we can get 


nQo(go,q,7',™o,r) = 


me 


2 rl 


dx 


_f]_ 

P Vo ““’7 (27r)2 

[(1 - 2x)qo + d] S'*} , 


nii(go,q,T,mo,r) = - 


8 iVe' 


2 /.i 


dx 


dH 

(2^)2 


Si-2 


Si- 


P + mg + x{l - x)ql + a; ( /cq + -go ) ^ 


1^ + 2a;(l - a;) ( go + ^ ) + a; (2fco + 3go) 5 


[(1 - 2x)go + d] S*} , 


with 


s. = E 


1 


S' = E 


—00 W + mg + x{l - x)q^ + 2x {k^ + go) S] 

^0 


—00 [^0 + + mg + x{l - x)q^ + 2x {k^ + go) S] 


2 ■ 


(A4) 

(AS) 

(A 6 ) 

(A7) 


where P = Iq + with Zq = + xqo + rsgn(fco), g^ = go + q^j and d = E [sgn (fco + go) — sgn(fco)]- When d ^ 0, the 

frequency summation cannot be carried out precisely. There are two ways to make d = 0. First, go = 0, corresponding 
to the static limit. Second, E = 0, corresponding to the clean limit of the system (without any disorder). Next we 
calculate the polarization functions in these two cases respectively. 


2. Calculation of polarization functions in the limit go = 0 

For a general constant E, we have 


4Ap^ f /■ d^l 

noo(q,T,mo,E) = J dx J [Si - 2 (l^ + m^) S 2 ] , 

oi\Tp2 pi p j 2 i 

nj,(q,r,mo,E) =- — J dx J [Si - [l^ + a;(l - a;)q^] S 2 ] 


with 


- It,) E 


n=o [^(n + 5 + Xsgn (n + \))^ + 


where X = ^E and Y = + ttiq + x(l — a;)q^. Summing over n, it is easy to get 


Si = 


27r2r 


Im 


^(-+X + iY 


(A 8 ) 

(A9) 


(AlO) 


(All) 
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which then leads to 


c /3" 5^1 /34 ^ 

^2 — —— Trr~: , — tt— ^^ Im 




P'^ dlm[p:{\ + X+ iY)\ 


87r2y dY 167r4y3 

Substituting the expressions of Si into Eqs. (IA8I) and (IA9I) . we have 


167r4y2 


dY 


2iV6^ /*^ f c‘^ 

noo(q,r, mo,r) = —^ dx dt < -^Fiit,T,r) + 

2Ne^ /■! ( 

nii(q, T, mo,r) =-5— dx _ dt < 

^ Jo J x/rni+Cj: I 


-Cl dF^{p,T,T) 




t dt j ’ 

— TOq SlniEi (t, T, r) 


t 


dt 


(A12) 

(A13) 

(A14) 


where Cq = xjx[l — x)ap and E'i(t, T, F) = Im \4>{\ + 2 ^ + * 2 ^ 7 )] ■ oxe interested in the limiting behavior 

of IIoo and 11^^ at zero temperature. As T —>-0, we know that 


lim Im 

T->0 


,1 r + it 

^^ 2 ^ 2ttT 


= arctan I p I ■ 


(A15) 


Therefore, at zero temperature the polarization functions can be written as 

A \ 


2 iVe2 

noo(q,mo,r) = —^<rin 






2|q| ViF 2 +|q| 


Jo 


a;(l — x) f Ki 

dx --arctan | 


^ 2iVe2 2 f , cell-a; /Ki 

nj,(q,mo,r) = —^q / dx ---arctan — 

Jo Ki V r 


with Ki = y^tuq + a;(I — a;)q2 and K 2 = \/4 (r2 + tuq) + q2. 


(A16) 


3. Calculation of polarization functions in the clean limit F = 0 

For general goi we have 

4iVe^ /•! f ffi\ 

^oo{qo,<l,T,mo) = J J {'S'l - 2 [ 1 ^+ mo+ a;(I - a;)go] 52 + (1 - 2a;)go5'*} , (A17) 


0 

8Ne'^ f 

Ilii{qo,q,T,mo) = - — J J 


dH 

(27r)2 


5i- 


1 + 2a:(l - x) ( + 


,2 , q 


2 \ 1 


with 


5. = 


2i 00 


E 




n— — oo 

3 cx) 


[(n+i+A)%y 2 

n + ^ + A 


n—— 00 


2 ■ 


52 + (I - 2x)go5* [>. (A18) 

(A19) 
(A20) 


(n + i+A) +y2 

where A = -^^xq^ and Y = ^a/12 + TOq + x(I — x) {q^ + q2). Carrying out the frequency summation yields 


/3" 


n \i 


r /I \1 

47r2A 

il) 


47r2A 



(A21) 


It is then straightforward to obtain 


52 = - 


S* = - 


dS, 

P^ 

Stt^Y dY 

327r4A2 

dlm[tlj{^ 

- A + iY)] 

dY 

P dSi 

i ' 

47r dX 

167r3A 


-Im 


i/; ( - + A + iA 


5Im [^(i + A + iF)] I 


dY 


Y 


Im 


__A + iA 


,03 j gini [^(i + X + lA)] 9Im - A + lA)] 


dX 


dX 


(A22) 

(A23) 
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Substituting Eqs. (IA21|) . (IA22I) . and (jA23l) into Eqs. (IA17P and (IA18I) . the polarization functions can be written as 


Ne'^ f 

noo(qo,q,'r,mo) = — 5 - dx 

^ Jo J ^ 


dt ■ 


0 dv^mg+CI 
(1 - 2x) qo d{F 2 {x, go, t,T)) 


^F2{x, qo, t, T) + — - - - 


2 

2 pi 


d{xqo) 


Ne f f 

Iluiqo,n,T,mo) = - 5 - dx 

Jo j. 


dtt ■ 


F - B'/ + m, 


0 


t 2 


[F2{x,qo,t,T)) + 


(A24) 

e + B';^-mld{F2{x,qo,t,T)) 


dt 


( 1 - 2 .), 


(A25) 


where Bg = \/x{l - x)q^, B'g = yjx{l - x)ql, and t = + ttiq + a;(l - x) {q^ + q^). Here, F 2 {x,qo,t,T) = 

Im + 12^)] + Im + * 2 ^)] ■ Since V^(l - z) = ipiz) + 7rcot(7rz), we have 


F2{x,qo,t,T) = 


2i 


— tan 


xqo _ . J_ 
2T ~ * 2r 


tan 


xqo t 
2T ^^2T 


Substituting Eq. (IA26I) into Eqs. (IA24I) and (IA25|) . the polarization functions become 


Ne^ r (2 [ 

noo(9o,q,7’,™o) = - 1 ^ J 


1 


cosh I -fJK^ ) — sin -Pxqo 


1 


1 [mg + a;(l - a;)qQ] sinh(/3A:3) 1 {1 - 2x)qoSm{/3xqo) 


K 3 cosn (i/SAla) — siir{^(3xqo) 2 cos\i {\l3K^) — sin^(i/3a;go) 


Bii{qo,(\,T,mo) = 


Ne- 


2 /-i 


dx ■ 


1 x(l - a:) (2gg + q2) sinh(/ 3 A: 3 ) (1 - 2x)go sin(/3a:q'o) 


27r Jo 1^3 cosh (i/JATs) - sim(i/3a;qo) cosh (i^SATs) - sim(i/3a;go) 


with K^, = \/ml + x(l - x) {q^ + q^). 

At zero temperature, the polarization functions are simplified to 


noo(go,q, mo) = 


nii((7o,q,mo) = 


iVe^q^ 


mo q^ — 4mo 

V 


4g3 


■ arcsm 


\/g^"+4^ 


Ne^ (2gg + q^) 


mo g2 _ 477 j2 

V 


4g3 


■ arcsm 


In the limit go = 0, the polarization functions are 

.2 ri 


Ne^ f f 

noo((l,T,mo) = —^j da;|2rin 


Ne"^ 

Il,i{q,T,mo) = - / 

TT Jo 


dx 


a;(l — a;)q^ 


2 cosh I 

2T 


A/g2 + 4mg 


mn , I Ki 

- 2 tanh — 


Ki 


2T 


tanh 


(A26) 


(A27) 

(A28) 


(A29) 

(A30) 


(A31) 
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